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Abstract: We present a real-time method to determine the beam prop- 
agation ratio of laser beams. The all-optical measurement of modal 
amplitudes yields parameters conform to the ISO standard method. The 
experimental technique is simple and fast, which allows to investigate laser 
beams under conditions inaccessible to other methods. 
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1. Introduction 

The characterization of laser beams is traditionally a basic task in optical science and performed 
for decades now. However, the existence of different archaic approaches like the variable aper- 
ture or moving knife-edge method, which are known to produce deviant results, has shown the 
need for a standardization of the definition. Especially the question how to reproducibly and 
reliably measure beam quality, which is very important from an user-oriented point of view, 
stimulated a lot of discussions HI. The ISO standard 1 1 146-1/2/3 has brought a great unifica- 
tion by defining all relevant quantities of laser beams including instructions how to perform the 
measurements conform to the ISO standard l2]. In the general case of so-called general astig- 
matic beams this approach relies on the determination of ten independent parameters, which 
are the second order moments of the Wigner distribution function. From these, three quantities 
can easily be derived according to ISO 1 1 146-2: the beam propagation ratio simply called 
parameter, the intrinsic astigmatism, and the twist parameter 121 . Especially the factor has 
become a well accepted measure for beam quality in the laser community by now — although it 
has to be used with care. 

Siegman pointed out that the use of measurements not conform to the ISO standard will 
result in values for M^, which are not comparable to each other HI. This is an important fact 
since other techniques as the aforementioned knife-edge method are still in use for one reason: 
their simplicity. As stated above, a measurement of a general astigmatic beam, which is fully 
conform to the ISO standard, requires the measurement of ten second order moments of the 
Wigner distribution function and, thus, is experimentally cumbersome and slow. 

Fortunately, due to their high symmetry most laser beams of practical interest require less 
than these ten independent parameters to be measured for a complete characterization: Most 
lasers emit beams that are simple astigmatic or even stigmatic because of their resonator de- 
sign. In this case the determination is based on a caustic measurement that includes the 
determination of beam width as a function of propagation distance. To do this, the beam width 
has to be determined at least at 10 positions along the propagation axis, completed by a hyper- 
bolic fit. According to ISO 1 1 146-1, the beam width determination has to be carried out using 
the three spatial second order moments of the intensity distribution 

Despite its experimental simplicity, the caustic measurement still is quite time-consuming 



and requires a careful treatment of background intensity, measuring area, and noise, which 
makes high demands on the temporal stability of a cw laser or the pulse-to-pulse stability of a 
pulsed source. Caustic measurements are therefore unsuitable to characterize the fast dynamics 
of a laser system. To react on the need for a faster and more detailed analysis, several other 
methods for laser beam characterization were presented such as Shack-Hartmann wavefront 
analysis lO, measurement of the complete 4D Wigner distribution function p^, or the use of 
diffraction gratings Q. 

In this paper, we investigate the possibility to monitor the beam quality based on a decom- 
position of the laser beam into its constituent eigenmodes. Many laser resonators possess rect- 
angular or circular symmetry. In this case, the laser beam emerging from the resonator can 
completely be described as a superposition of the well-known Hermite-Gaussian or Laguerre- 
Gaussian modes From the point of view of beam quality, they represent a natural choice 
of description since the fundamental Gaussian mode has the lowest possible value of = 1 . 
Any deviation from the ideal diffraction-limited Gaussian beam profile can be attributed to the 
contribution of higher order modes, leading to > 1 . Here, any excited higher order mode 
contributes with a certain value to the beam propagation parameter. This value can be readily 
calculated for every mode in a way which is conform to the ISO standard 13. Therefore, per- 
forming a modal decomposition of a given beam enables us to determine a value for which 
is compatible to the one resulting from a caustic measurement. 

To determine the modal weights necessary for the decomposition, various methods have been 
suggested such as coherence measurements |8 |, intensity recordings [9-11], or the use of a ring 
resonator |T2l. Again, one has to see the required experimental or numerical effort that makes 
those methods unsuitable for the intended task of real-time characterization. Hence, we suggest 
to use the correlation filter technique based on computer-generated holograms (CGHs) ifTslfT?! . 
Recently we successfully applied this method to characterize the field produced by multimode 
optical fibers and to determine the polarization states of different modes 10311161 . 

Here we will show that the correlation filter method yields values conform to the ISO 
standard but with diverse advantages. The experimental realization is simple and yields in- 
formation about the beam quality in real-time in contrast to other approaches. It relies on a 
decomposition of the electromagnetic field into the spatial modes of the resonator generating 
the laser beam. Since our investigated laser cavity possesses rectangular symmetry, Hermite- 
Gaussian modes are used in this work, but represent no restriction of the method. Any complete 
set of suitable modes may be used in other cases. 

From the measured modal weights, we calculate the beam propagation factor and com- 
pare our results to values obtained from traditional caustic measurements as defined by the ISO 
standard. Moreover, we demonstrate the on-line monitoring of the modal spectrum as well as 
the factor while realigning the resonator in real-time. 

2. Modal decomposition 

2.1. Expansion into Hermite-Gaussian modes 

Laser resonators with rectangular geometry generate superpositions of nearly Hermite- 
Gaussian (HG) modes with the field distribution ifTTl 

HG,„„(x,y) = -\/ ,„,+„^ , , H,„ f ^x") H„ \ , (i) 

Wo V 2"'+" TTOTln! \ Wo / Y Wo / \ J 

where wo is the waist radius of the fundamental mode HGoo and H; denotes the Hermite poly- 
nomial of order /. 

Hermite-Gaussian modes are one set of orthogonal eigenfunctions of the scalar Helmholtz 



wave equation. Due to the completeness of this eigenfunction set, an arbitrary transverse wave 
field U can be expanded into a superposition of HG modes ||6l 



U{x,y) ^Y.Y. c,„„HG,„„(x,y) , c,„„ = // \lGl„UAxdy , (2) 

m=On=0 JJ-^ 

where the asterisk denotes complex conjugation. The complex-valued expansion coefficients 
Cmn = Pmn^"^'"" include modal amplitudes p,„„ = \c,„„\ and phases ^„,„ = arg(c„j„). Using nor- 
malized fields with unit power P = Jf^^ \U\^ dxdy = 1, the squared amplitudes p,^„, represent 
the relative power of the HG,„„ mode since „ p,^,,, = 1 . 

2.2. Beam propagation ratio 

According to the ISO standard Q, the beam propagation ratio for simple astigmatic beams is 
defined by 

where d is the beam diameter and denotes the divergence angle. These parameters are based 
on second order moments C7^. For instance, the beam diameter in x-direction reads as i2j| 



Using HG modes, the beam quality from Eq. ^ simplifies to 
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In other words, the ratio of the beam waist diameter to the one of the fundamental Gaussian 
beam determines the factor. 

The laser used for experimental demonstration (manufactured by Smart Laser Systems) gen- 
erates simple astigmatic beams composed of different HG,„o modes with higher order modes 
in one direction ifTSl . In general, different modes of a resonator possess slightly dissonant fre- 
quencies ifTTll . Only modes with p = m + n ^ const belonging to one mode group have the 
same frequency and, thus, temporally stable intermodal phase differences. Due to the huge 
camera integration time compared to the beating period, the interference terms between modes 
of different mode groups vanish, whereas modes of one and the same mode group contribute 
coherently to a recorded intensity distribution. Furthermore, a rotation of the principal axis of 
one mode towards the laboratory system can be described by a superposition of HG modes 
of the associated mode group — similar to the expansion of Laguerre-Gaussian modes into HG 
modes ||6|. Consequently, the total intensity distribution is the sum of intensity contributions 
from single mode groups 

p=0 k,m 

where p denotes the /?* mode group and Ip is the intensity formed by a superposition of the 
modes from mode group p. 

Using the intensity profile from Eq. (|6]l, the first order moments of the beam vanish and, thus, 
express that the beam centroid is on the optical axis. Hence, the spatial second order moments 
read as ^ 

-2- fl'°°^x^I{x,y)dxdy=^Y.pl,{2m+l) , (7) 
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and analogously for a^, where the identities lfT9l 
HG|(./HG,„„ djcdy : 
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xyHGwHG,„„dxdy = -^^max(^,m) 5|i._,„|,i Vmax(Z7;) 5|/_„|_i , (10) 

and the Kronecker symbol 5,/ have been used. 

Combining Eq. (|4]i and (|5]l with Eq. ^ and (HJl, the beam propagation ratio in jc-direction 
reads as 



Ml = ^pi(m + « + l) 



r 



\m.n J 



+ 1 + i 



1/2 



(11) 



To obtain , the plus sign in front of 7 has to be replaced by a minus. Due to Eq. (fTTl i. the 
knowledge of the modal weights is sufficient to determine the beam propagation ratio . 

3. Experiments 

3.1. Measurement setup 

The optical setup in Fig.[T]consists of a laser source and two branches that enable two simulta- 
neous but independent measurements of the parameter of the emerging beam. The Nd: YAG 
laser can produce different HG„,„ modes at A = 1064nm. The excited HG,„„ modes are re- 
stricted to n = since the pump light (A = 808 nm) is coupled into the plane-concave resonator 
by a horizontally movable fiber ifTSl . A beam splitter divides the beam and guides the replicas 
into the two branches for analysis. 
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Fig. 1. Scheme of experimental setup. Upper branch: M determination using a CGH. 
Lower branch: determination by a caustic measurement conform to the ISO standard. 



In the first branch a lens images the laser beam waist onto an adapted CGH, which is an 
amplitude hologram consisting of 601 x 601 Lee cells with cell widths of 16 \im. The far field 



intensity behind the CGH realized by a subsequent 2/-setup is recorded by a CCD camera. This 
intensity pattern contains the information about the weights p,^,„ of the HG modes since the 
complex conjugated modes are implemented in the CGH called correlation filter jTSl . Hence, 
the correlation filter technique enables the determination of the factor in real-time using a 
computer-aided calculation of Eq. ( fTTT i. 

The CGH for this experiment is designed to simultaneously analyze the amplitudes of 21 HG 
modes, i.e., all HG„,„ modes with (m + «) < 5 are implemented taking into account a possible 
rotation of the resonator coordinate system. This correlation filter configuration corresponds to 
a truncation of all sums in the equations of section |2] 

In the second branch of the setup in Fig.[T] an additional lens is used to analyze the caustic 
of the beam conform to the ISO standard 121, which serves as a reference measurement. 

3.2. Measurement results 

To investigate the reliability of the CGH-based measurement procedure, a modal spectrum of 
an arbitrarily chosen beam consisting mainly of HG20 and HG30 is depicted in Fig. |2] Using 
Eq. (|6]l the intensity distribution in the plane of the CGH can be reconstructed. Since the method 
will only be as good as it is able to reconstruct the investigated beam, the measured intensity 
distribution is compared to the reconstructed one by calculating their two-dimensional cross- 
correlation coefficient. A value of 1.0 denotes perfect match. Cross-correlation coefficients 
greater than 0.9 in all investigated cases indicate the excellent functionality of the method. 
In particular, the cross-correlation coefficient for the beam investigated in Fig.|2]is 0.98. For the 
determination of M^, the reconstruction of the intensity profile is not necessary as a matter of 
course and just shown here to illustrate the functionality of the method. 
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Fig. 2. Modal spectrum of an arbitrarily chosen mode mixture (left) with corresponding 
measured and reconstructed near field intensity distributions (right). 

The results for the measurement of the beam propagation ratio using the correlation filter 
technique and the ISO standard method are compared in Fig. [3] Seven different mode mixtures 
are selected and analyzed to test the accordance of the results of both techniques for beams with 
different factors. For instance, the beam in Fig. |2]is denoted as Mix 6. As the bar charts in 
Fig. [3] show, the results of both measurement techniques are in very good agreement. The M^r 
parameters confirm the characteristic of the laser to emit only mixtures of HG„,o modes consis- 
tent with the theory. The slight deviations in the factors are caused by several reasons: First, 
the hologram needs to be placed accurately at the waist position and on the optical axis. Fur- 
thermore, the waist diameter has to match the beam diameter the correlation filter is designed 



for. Any deviance from the ideal alignment as well as CCD noise affect the measured modal 
weights resulting in slight inaccuracies of the CGH-based beam quahty determination. 
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Fig. 3. Comparison of the Mj and My factors determined by the CGH technique and the 
ISO standard method, respectively. 

The deviation between the differently determined factors behaves statistically, whereas a 
tendency in the deviation of the values is recognizable: For mixtures with low factors 
such as Mix 1 to Mix 4, the values determined with the CGH are larger than the ones of the 
ISO conformable measurement. Here slight errors of the modal weights occur for higher order 
modes due to CCD noise producing an offset of the implemented higher order mode amplitudes, 
which enlarge the CGH-based factors. On the other hand. Mix 6 and 7 mainly consist of 
modes from mode groups with m + n near the truncation limit. Due to the finite number of 
implemented modes in the CGH, modal weights of HG,„„ modes with (m + n) > 5, which would 
increase the beam quality Mj:, are set to zero in the experiment. Hence, the missing amplitudes 
cause reduced factors. For MixS these two effects nearly compensate one another. 

In all analyzed cases, the maximum measured deviation is 13% forM^ and 5% forM^. This 
result demonstrates the functionality of the method, especially when beams with high beam 
quality are considered. 
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Fig. 4. mI and My on a time scale of 30 s taken with a rate of 4 Hz. To vary the beam 
quality the laser resonator was tuned continuously. The insets depict the waist intensities at 
selected points in time. 



To demonstrate the real-time capability of the all-optical CGH technique, we continuously 
varied the alignment of the laser cavity and monitored the output. Figure |4] illustrates an hfi 



measurement on a time scale of 30 s, where the cavity was initially aligned to produce the fun- 
damental HGoo mode. Then we changed the transverse position of the end-pumping fiber ifTSl 
and the factor increased. After 15 s we reversed the process to return to the fundamen- 
tal Gaussian mode after 30 s. The reachable speed of the method is only limited by the used 
hardware. 

4. Conclusion 

We have shown that a modal decomposition of an investigated laser beam allows not only the 
quantitative determination of its transverse modal content but also of its beam quality factor 
M^. Using the correlation filter technique, it is possible to determine the modal spectrum in 
real-time. We demonstrated the working capabilities by analyzing the beam emerging from an 
end-pumped Nd: YAG laser By changing the alignment of its cavity, several Hermite-Gaussian 
modes of higher order could be excited. We showed that a correlation analysis involving the 
first 21 Hermite-Gaussian modes yields values for a derived M" parameter, which are in good 
agreement with the caustic measurements proposed by the ISO standard. The advantage of the 
method is the ability for real-time analysis of the beam that was demonstrated by continuously 
changing the alignment of the laser cavity and monitoring the beam propagation ratio as a 
function of time. 

We conclude that this method provides an experimentally simple possibility to measure the 
beam propagation ratio M^, which is compatible to the ISO standard in cases where caustic 
measurements are too time-consuming or even impossible to perform. The additional informa- 
tion about the modal spectrum of the beam is highly advantageous for diverse laser applications. 



